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Abstract. The J = 2 — 1 and 3 — 2 rotational lines of the rare isotopomer C^^S and the J = 5 — 4 and 6 — 5 
transitions of Si^'^S were detected in the carbon star IRC+10216 (CWLeo). These are the first detections of S 
bearing molecules in a star and the first spectroscopic detection of Si"^^S. From a comparison of ■^'*S and "^^S 
bearing isotopomers, the ^*S/'^''S isotopic ratio is 107(±15). This value is comparable to values in the interstellar 
medium of the inner Galactic disk (115) but is smaller than the solar value of 288 (Ding et al. I2001II . The 
increase of the ^®S abundance relative to ^''S only qualitatively follows model predictions of a low mass AGE 
star. Quantitative agreement of the observed "^^S/'^^S ratio with the stellar models can be reached if the age 
of IRC+10216 and Galactic chemical evolution are taken into account. Other less likely possibilities are the 
presence of considerable inhomogeneities in the interstellar medium and either IRC+10216 or the Sun started 
with a peculiar '^^S abundance. Other production mechanisms potentially capable of enhancing the Galactic 
interstellar medium are discussed. From the observed line density toward IRC+10216 and toward Galactic star 
forming regions, we estimate the confusion limit toward those sources. 
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1. Introduction 

Heavy elements are formed in stars of medium or high 
mass. At the end of their lifetimes, these stars recycle, via 
winds or via explosions, a fraction of these metals back 
into the interstellar medium (ISM). Over the aeons, inter- 
stellar abundances of such heavy elements have been in- 
creasing and isotopic compositions have been changing, as 
can be seen by comparing the composition of the present 
day ISM with that of Damped Lyman-a systems (Rauch 
I1998|l . Not only elemental but also isotopic abundances, 
which are mainly determined in the radio range, are pow- 
erful tools to investigate stellar production sites and the 
star forming history of galaxies (Wilson & Matteucci ll992l 
Henkel & Mauersberger Henkel et al.raM!, Wilson 

& Rood 1994, Kahane 1995, Prantzos et al. il996i Bieging 
[TM7I Henkel et aL lTMSji . 
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Being one of the ten most abundant elements and pos- 
sessing four stable isotopes, sulfur is of particular inter- 
est for such studies. Relative abundances of stable sul- 
fur isotopes have been determined for meteorites (e.g. 
Gao & ThiemenslHnil), the Moon ( Thode & Rees lWTIl . 
the Galactic ISM (e.g. Chin et al. HM^ . Cosmic Rays 
(Thayer I1997|l and late type stars (e.g. Kahane et al. 
11988(1 . Even in external galaxies ^^S/'^^S ratios could be 
measured (Mauersberger & Henkel [T989I Johansson et al. 
I1994f) . Small isotopic variations of sulfur in terrestrial, me- 
teoric and planetary samples have been proposed to be 
powerful indicators of chemical, geophysical and biologi- 
cal processes (Canfield I2UU1I Farquhar & Wing 2003 ) 

Mauersberger et al. H199(i(l presented the first interstel- 
lar detections of a ^^S bearing molecule toward a num- 
ber of Galactic hot cores. While the solar '^^S/'^^S ratio 
is 288, an interstellar ratio of 115(±17) was found. This 
supported the notion that unhke other S isotopes, ^^S is a 
secondary-like nucleus, predominantly synthesized via the 
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Fig. 1. The whole spectral range measured near the C^^S 
J = 2 — 1 (lower panel) and 3 — 2 (upper panel) transitions 
toward IRC+10216. The channel spacing corresponds to 
3.1 kms"^ for the 3mm spectrum and 4.2kms~^ for the 
2 mm spectrum. The intensity is in mK (T^). For identi- 
fications, see the AppendixIXI 

s-process in massive stars. Most interstellar '^^S/'^^S ratios 
were determined, however, for sources with galactocentric 
distances Rqc < 7kpc. A positive '^''S/'^^S gradient with 
increasing RqCi i-e. a ratio between 115 and ~ 180 in the 
solar neighborhood, is not completely ruled out by the 
data of Mauersberger et al. (1996). 

In order to further constrain the production site of '^^S 
in the Galaxy, we conducted a search for rotational lines 
of C^^S and Si'^^S toward the prototypical carbon star 
IRC+10216, where all other stable S isotopes have been 
detected previously (Kahane et al. I1988|l . 

2. Line frequencies 

The rest frequencies of the C"^^S J = 2 — 1 and 3 — 2 lines 
are 95.0167 and 142.5228 GHz (V. Ahrens, pers. comm.), 
and 96.4130 and 144.6171 GHz for the corresponding C^^^S 
lines (Lovas 1992). 

The frequencies of Si^^S have not been measured in 
the laboratory. However, it is straightforward to derive the 



rotational constants of Si'^^S from those of SiS (see, e.g., 
Frum et al. HMIjl using the isotopic relation for diatomic 
molecules (e.g. Townes & Schawlow I1975|l . For Si^^S, 
we derived Bo=8607.47MHz and Do=5.4 lO'^ MHz. The 
computed frequencies are i^(J = 5 — 4)=86071.984MHz 
and t/(J = 6 - 5) = 103284.96MHz. The expected relative 
error for the rotational constants is < 10^^ and the re- 
sulting expected uncertainty for the J = 5 — 4 and 6 — 5 
lines of Si'^^S is < 0.5 MHz. We have checked the preci- 
sion of those relations for the different isotopes of SiS that 
have been measured in the laboratory and we find that 
the quoted uncertainties are rather conservative. Recently, 
Sanz et al. |2003| measured the rotational constants of 
several isotopes of SiS and have provided a fit, taking into 
account the breakdown of the Born-Oppenheimer approx- 
imation, to the data of all isotopes of SiS. From their con- 
stants we infer those of Si^^S to be Bo=8607.495MHz, 
1)0=5.364 10-3 MHz, jy(J = 5 - 4)= 86072.27MHz and 
jy(J = 6 - 5) = 103285.31MHz, i.e, 0.3 MHz of difference 
with respect to our early calculations. Hence, we are con- 
fident that the calculated frequencies have an error well 
below 1 MHz. 



3. Observations 

The J = 2 - 1 and J = 3 - 2 lines of C^^S and C^^S were 
observed between October 1999 and September 2003 with 
the IRAM 30 m radiotelescope on Pico Veleta (Southern 
Spain) toward IRC+10216. The Si^'^S observations were 
carried out between 1996 and 2003 at the same telescope. 
Several spectra at different frequency settings were ob- 
served for each line of Si^^S to ensure that the observed 
lines are coming from the signal side band. The observa- 
tions were conducted under good weather conditions. SIS 
receivers were used with image sideband rejections of 20 — 
25 dB. An antenna temperature (TJ) scale was established 
by a chopper wheel method. The forward efficiency of the 
telescope was 0.94 and the beam efficiencies were 0.80 at 
3 mm and 0.65 at 2 mm. The beamwidth (full width at 
half maximum) of the 30 m antenna was 26" and 17" for 
the 2 — 1 and 3 — 2 transitions of the isotopomers of CS; 
for the 5 — 4 and 6 — 5 transitions of SiS, beamwidths 
were 28" and 24". From scans of continuum sources, we 
estimate the pointing to be correct within 5". 

As backends, we employed filterbanks with 256 chan- 
nels and a channel spacing of 1 MHz corresponding to 
~ 3.1 kms"^ for the 3mm transitions and 2.1 kms^^ for 
the 2 mm transitions. Two receivers were employed simul- 
taneously, either observing orthogonal polarizations at the 
same wavelength or measuring one line at 2 mm and the 
other at 3 mm using just one polarization for each wave- 
length. 

All spectral lines were measured using the wobbling 
secondary mirror with a beam throw of 200" in azimuth. 
The phase time was 2 seconds and the on-source integra- 
tion of each subscan was 30 seconds. Subscans obtained 
by wobbling to the left and to the right of our source were 
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Table 1. Line parameters 
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a) the ratio between the intensities at the central velocity of 
the line and at the velocity of the horns (i.e. maximum blue or 
red shifted emission); b) fixed to values obtained for the C^^S 
or Si'^'^S lines; c) simultaneous fit with a blended, unidentified 
line at u = 103.2804GHz and an intensity of 46(18) Kkms'^; 
d) from the integrated intensities yielded by the fits to the lines 
of ^*S and ^®S bearing species. The error given is the 1 a error 
from the fit and does not contain an estimate of systematic 
errors, e.g. due to baseline uncertainties. 

added to eliminate baseline ripples caused by the asym- 
metry in the beam path (symmetric switching). 

4. Results 

Fig. n displays the entire spectral range measured at the 
frequencies of the C^^S Hues. Figs .Hand El show the C^^S, 
C^^S, Si^^S and Si^^S spectra that are discussed in Sect. 
15.21 in more detail. An inspection of Fig.^ shows that a 
main source of uncertainty in the identification of weak 
lines and also in the definition of spectral baselines comes 
from the presence of many weak line features, in particular 
at 2 mm. In AopendixIXl we list the line parameters of all 
features together with possible identifications. The Si'^^S 
data, showing a statistical behavior that is similar to the 
C^^S data, will be analyzed elsewhere (Cernicharo et al., 
in prep.). 

From the spectra, we subtracted baselines of first or- 
der. For the 3 — 2 line of C'^^S, almost the entire frequency 
range is covered with line features near or above the de- 
tection threshold making the identification of a region free 
of emission a difficult task. We chose spectral regions with 
particularly low emission for the definition of a baseline. 

5. Discussion 

5.1. Where is the confusion limit toward IRC+10216? 

As it can be seen in Fig.Q] showing the total spectral 
range of 2x256 MHz toward IRC-H10216, we have de- 



tected a large number of spectral line features, many of 
them unidentified. There are more lines at 2 mm than at 
3 mm although the rms of the data is similar. The spec- 
tral line density in our 2 mm spectrum is 15 lines per 
256 MHz. Since each line has a width of about 14.5 MHz, 
about 85% of the spectral range observed would be cov- 
ered with detectable line emission if there were no overlap. 
The weakest feature we identified in Table TOI at A=2 mm 
has r^=3.5 mK, although the rms noise of our 2 mm spec- 
trum is about 0.5 mK for 4 kms~^ wide velocity channels. 
Apparently lines weaker than about 3.5 mK have escaped 
detection not because of lacking signal-to-noise but be- 
cause their number is so high that they are all blended 
and therefore cannot be easily distinguished. 

At A = 3 mm, 6 lines have been detected from the sig- 
nal band in a 256 MHz wide range. Since each line has 
a width of 9 MHz, 21% of the spectral range observed is 
covered by detected lines. The weakest of those have a 
of about 1.5 mK, which corresponds to the 3a noise level 
at a velocity resolution of 10 km s~^ . In contrast to the sit- 
uation at 2 mm, we have not yet reached the 3 mm confu- 
sion limit. Nevertheless, at 3 mm we also see signs of weak 
blended lines at a level of ~lmK. Weaker lines, on the 
order of 0.5 mK, have been identified recently (Cernicharo 
et al. 12004(1 making use of the characteristic line shapes of 
IRC-H10216. 

To summarize, at the IRAM 30-m telescope the confu- 
sion limit toward IRC-l-10216 is reached at A=2mm after 
integrating down to a system noise of about 1 mK, while 
at 3 mm one has to integrate down to 0.3 mK. For other 
well known molecular sources, confusion limits are given 
in AppendixiBl 

5.2. The identification of (?^S and Si^^S 

As already mentioned in Sect.0]there are emission features 
at the rest frequencies of both the 3—2 and the 2—1 lines of 
C'^^S. These clearly peak above the limits of line confusion. 
It turns out, however, that at 3 mm there are two emission 
lines of linear C3H from the image sideband, one of which 
coincides with the location of the expected C^^S 2 — 1 
feature. Also in the band is the image of the strong CS 
2—1 line. From a well calibrated spectrum of these three 
lines (Mauersberger et al. I1H8 9»1 we obtained their relative 
intensities. According to Cernicharo et al. H2000|l there 
is no hint for variability in time, except for some very 
highly excited lines or maser lines. In our spectrum of 
C36s 2 - 1 (Fig.EJ we could therefore apply a fit to the 
CS line and the two C3H lines from the image band fixing 
the line shape and the relative intensities. When removing 
this composite fit from the spectrum, the residual still 
shows emission at the expected frequency of C^^S (Fig.|21). 
This residual spectrum was used to determine the isotopic 
ratios given below. 

We cannot exclude that the residual line is contami- 
nated since a fit with fixed line shape and rest frequency 
to the line feature (to the values obtained for C'^^S) shows 
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some excess emission toward lower velocities. A possible 
contaminant of the 2 — 1 line of C'^^S is a series of K 
lines of the J = 23 - 22 transition of CHgCCCN (methyl 
cyanoacetylene; Lovas IH^Mjl . which were detected toward 
a number of Galactic hot cores (Mauersberger et al. .l996 ) . 
This is, however, unlikely, since this molecule is not yet 
identified toward IRC+10216 (e.g. Cernicharoet aL lMjUjl . 

We have used the procedure from the CLASS data 
reduction package to fit shell type circumstellar lines. In 
order to obtain more reliable values for the integrated in- 
tensities, we fixed some of the line parameters of '^^S bear- 
ing species to the values obtained from the more intense 
^■^S bearing isotopomers. Fit results are given in Tablc^ 
In the fit to the C'^^S 3 — 2 line feature there is a hint for 
excess emission toward lower velocities. At that frequency, 
the only previously observed interstellar line in the Lovas 
HI 984(1 catalog is a transition of CH2CHCN (vinyl cyanide, 
V = 142.523 GHz). This complex molecule has never been 
detected before in IRC+10216 and the line strength is 
very low. It seems plausible that if there is contamination 
it comes from a transition of a molecule abundant in late 
type carbon stars but not in the ISM, such as C7II and/or 
H2C6 (Guelin et al. lTMTIl . 

The Si^^S 5-4 line can be weU fitted, with the de- 
termination of the baseline being the major source of un- 
certainty. The Si^^S 6 — 5 line is also clearly detected, 
although it seems to be blended with a weak, unidenti- 
fied feature (Fig.OJ. The parameters of the C^^S, C^^S, 
Si^^S, and Si^^S, lines and the derived C^^S/C^es and 
Si'^^S/Si'^^S intensity ratios are given in Tabled We can 
summarize that we have clearly detected C'^^S and Si'^^S in 
IRC+10216. These are the first detections of ^^S bearing 
molecules in a stellar atmosphere and the first detection 
of Si'^^S outside the solar system. 

5.3. The ^^Sf'^S abundance ratio 

The characteristic line profiles of spherically expanding 
molecular shells can be explained in terms of line opac- 
ity, expansion velocity and the relative sizes of the tele- 
scope beam and the stellar envelope (e.g. Kuiper et al. 
[BTfni Olofsson et al. [MS)- The shghtly U-shaped pro- 
files seen in C'^^S and Si'^'^S can be explained by opti- 
cally thin and spatially resolved emission. Optically thin 
lines are expected because the main isotopic lines are only 
moderately optically thick and because '^"'S is a factor 
~ 20 less abundant than the main isotope (Kahane et 
al. 119881 Cernicharo et al. 12000(1 . Interferometric observa- 
tions (Guelin et al. 119931 Lucas et al. 1995, for maps see 
Grewing 11994(1 confirm a compact condensation of about 
10" diameter for main isotopic CS. SiS is even more com- 
pact (Lucas et al. il995ii . which is also apparent by a com- 
parison of C^^S and Si^^S fine shapes (Figs. El and EI). It 
is safe to assume that also C'^^S and Si^^S emission is op- 
tically thin and it is plausible that the distributions of '^''S 
and '^^S bearing isotopomers are similar. Chemical frac- 
tionation, as discussed in Chin et al. ((1996(1 . should be 



negligible in the warm (50 K, Cernicharo et al. 12000(1 en- 
vironment of the inner shell of IRC+10216. Since the iso- 
topomers have very similar transition frequencies we can 
assume that the intensity ratios of corresponding transi- 
tions refiect the abundance ratios (for a discussion, see 
Kahane et aLgSHHIl- 

From Tabled it is evident that the relative errors of 
the integrated line intensities of the "^^S bearing molecules 
are of the order of 20% or more, while the relative errors 
for the ■^'^S bearing species are much smaller. While it 
is relatively straightforward to determine the X'^^S/X^'*S 
ratios from fitted line profiles (as summarized in Table^, 
it is much more difficult to estimate the X^^^S/X^^S ratios 
since the simple rules for error propagation apply only if 
the relative errors are small. In addition, if we assume that 
the error distribution of the X'^^S/X'^^'S ratios is Gaussian, 
the distribution of the inverse ratio is not such a simple 
and symmetric function, making the determination of the 
variance difficult (see Chapter 3.3.3 in Wall & Jenkins 
I20Q3.) . 

The errors as quoted have been calculated from the 
formal errors of the fits to the line profiles. From these 
four ratios, and using as weights the inverse squares of 
the formal errors of the individual ratios, we have deter- 
mined the weighted mean of the '^^S/'^'^S abundance ratio 
and its error to be 0.0094(±0.0013) (see e.g. Bevington & 
Robinson I 1992(1 . From a graphical representation of the in- 
dividual values and their mean (see Fig.EJ the only transi- 
tion whose abundance ratio deviates from the mean value 
is the SiS 5 — 4 line. This might be an additional hint that 
for this line there is, in addition to the statistic uncer- 
tainty, also a systematic effect, for example a blend with 
an unidentified line (see the discussion above). Since the 
relative error of the average '^^S/^^S abundance ratio is 
only 14%, the expectation value of the •^"'S/'^^S ratio and 
its error can be easily determined. In the following we will 
use a 34S/36S ratio of 107(±15) for IRC+10216 to ease a 
comparison with values published in the literature. 

This value of the best fit is very similar to the inter- 
stellar value of 115(±17) (Mauersberger et al. 11996(1 but 
is smaller than the solar system value of 288 (Ding et al. 
12001(1 . In Table El '^e give a compilation of all available 
estimates for sulfur isotope ratios. 

5.4. The origin of^^S 

5.4.1. Solar abundance and nucleosynthesis of '^^S 

■^^S is one of the more enigmatic stable nuclei in na- 
ture. Partly because of its rarity (solar system isotopic 
abundance 0.015% relative to '^^S), and partly because it 
is rarely employed in mass spectrometric studies of sul- 
fur isotopes, up to very recently even its solar system 
abundance was quite uncertain. The often used Anders 
& Grevesse p989(l abundance table lists the S isotopes 
with a '^''S/^^S ratio of 211 (no errors given). An up- 
dated value is recommended in the 1998 compilation of 
isotopic abundances of the International Union of Pure 
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Table 2. Isotopic ratios of sulfur" 



Ratio 


solar 


ISM 


C.R. 


IRC+10216 




system 










22.64(0.00) 


24.4(5)" 


16(-5,+ll) 


21.8(2.6) 


32g^33g 


126.95(0.05) 


~150 


38(-18,+500) 


121(15) 


32g/36g 


6519(20) 


3280(760) 




2700(600) 


34g/33g 


5.606(0.002) 


6.3(1) 




5.6(.3) 


34g/36g 


288(1) 


115(17) 




107(15) 



a) References: solar system data; Canyon Diablo troilite (CDT), Ding et al. 120011 ISM data: Chin et al. (I1996II . Mauersberger et 
al. ( 1996), most ISM data are from the inner Galaxy; Cosmic Rays: these are "source abundances", Thayer 119971 : IRC+10216: 
Kahane et al. 1200011 . this paper. 

b) the average value for sources at various galactocentric radii; the data by Chin et al. 119961 indicate a possible galactocentric 
32gy34g abundance gradient, which would correspond to a value of 32(±5) at the solar circle. 



and Applied Chemistry, (Rosman & Tavlor n"998|l , where a 
ratio ^"'S/^^S = 248±29 is hsted as the "best measurement 
from a single terrestrial source" ; this composition has also 
been used by Lodders H2003|l for her Solar System abun- 
dance table. Sulfur isotope ratios measured in the labora- 
tory are commonly reported as permill deviations from the 
Canyon Diablo Troilite (V-CDT) standard, which has an 
S isotopic composition fairly representative for Solar sys- 
tem materials (TUPAC 2003 compilation; de Laeter et al., 
I2()()3|l . Its composition has recently been precisely redeter- 
mined, with a 34S/36S ratio of 288±1 (Ding et al.,||Mil|)- 
Another modification with respect to Anders & Grevesse 
l|1989|l comes from a lower sulfur elemental abundance by 
a factor 0.86, as first deduced by Palme & Beer (1993) and 
quoted in Lodders (|2I)03) for the solar system. All in all, 
the solar abundance of ^^S is lower than the commonly 
used Anders & Grevesse (^12S2J value by a factor 1.6. 

The former uncertainty on this isotope extends to its 
nucleosynthetic origin. The synthesis of the major isotopes 
of sulfur, similar as for other elements between Si and Fe, 
can mostly be assigned to charged particle reactions dur- 
ing hydrostatic and explosive burning stages in massive 
stars: '^^S and '^'^S to hydrostatic and explosive oxygen 
burning, ^^S to explosive oxygen and neon burning (Chin 
et al. 1996; Woosley et al. 2002). For some exceptionally 
rare neutron-rich isotopes in this mass range including 
^^S, however, the situation may be different. ^^S has a 
closed neutron shell, hence a small neutron capture cross 
section, and in an s-process it could, in principle, act as a 
bottleneck building up a significant abundance relative to 
neighboring nuclides. For this reason, the weak s-process 
taking place during core He burning and convective shell 
C (Ne) burning in massive stars has been suggested as a 
major source of the ^^S abundance in the Solar System 
(Schatz et al. 119951 Woosley et al. l2002|l . and according to 
first network calculations by Schatz et al. H1995(l it seemed 
to be able to fulfill that role. However, re-assessment by 
Reifarth et al. H2000|l with updated nuclear input data in- 
dicates a more complex situation. The unexpectedly low 
neutron capture cross section of '^^S causes this nuclide 
to restrict the flow towards ^^S. With rare '^^Cl being the 
main seed nucleus rather than abundant ^*Si and ■^^S, the 



weak s-process, while actually making ■^^S, falls short by 
more than an order of magnitude in quantity in their as- 
sessment. Here, a further complication arises because in 
all recent calculations the abundance of the major seed, 
•^^Cl, which in the s-process zone is given by its original 
abundance, has been taken from Table 3 of Anders and 
Grevesse II198 9'I. Actually, there is disagreement for CI 
between the nuclide abundances in their Table 3 and the 
39% higher elemental abundance given in their Table 1 ; 
the latter agrees with other previous and more recent com- 
pilations. The predicted •^^Cl/'^^Cl ratio is not affected by 
these changes, but there is a significant effect on the ■^^S 
overproduction factor, which is increased by a factor of 
1.4. 

According to current models, in massive stars most 
■^^S is produced in a large mass zone where convective 
shell C-burning operates in hydrostatic conditions before 
the SNII explosion (Woosley & Weaver [l995|l . Owing to 
the high neutron density in the first phase of carbon con- 
sumption (Raiteri et al. I1992|l . more "^^S is fed through 
the neutron channel of the unstable '^^S than would be 
inferred by the occurence of the classical weak s-process. 
The nucleosynthesis yields of the most recent massive star 
calculations, made with solar initial composition and with 
a full network included (Woosley et al.|2002; Rauscher et 
al. nna^ Llmongl & Chieffi l^nMjl . are aU based on the 
Bao et al. H2()()()|l compilation of neutron capture cross 
sections (where the new and much lower cross section for 
■^■^S is included). These results show a production factor 
of "^^S that is typically low by a factor of 2 with respect 
to the mean of the most abundant nuclides. Taken at face 
value, this factor of 2 may be reconciled using the new so- 
lar ■^^S abundance as well as the correct '^^Cl abundance. 
However, this does not completely solve the problem of 
the reproduction of solar '^^S, since one would expect the 
weak s-process products to decline with decreasing metal- 
licity, i.e. earlier generations of massive stars to have been 
less effective producers of "^^S. In addition, there are un- 
certainties related to the nuclear input data. In particu- 
lar, the rate of the major channel '^^Cl(n,p)'^^S may be 
subject to further refinements in the future owing to the 
much better energy resolution obtained in the new mea- 
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surements of Wagemans et al. (|2Ul)3|l . On the other hand, 
recent investigations of n-induced reactions on "^^Ar and 
39 Ar (Goeminne et al. IMM I^HU^ seem not to affect the 
■^^S production in a significant manner. 

A complementary, though probably less important, 
source to the weak s-process is the main s-process tak- 
ing place in the He burning shell of low mass AGB stars. 
According to calculations partially published in Travaglio 
et al. H2()()4(l . the Galactic contribution to ^^S at the epoch 
of Solar System formation by all previous generations of 
AGB stars is about 4%; this may reach ^10% once the 
revised lower ^^S and higher '^^Cl solar abundances are 
taken into account. 

5.4.2. Stellar mass of IRC+10216 

In Fig.|Slwe show the S isotopic data for IRC+10216 (from 
this work and Kahane et al. I2000|l normalized to '^''S and 
solar system ratios. The data clearly show a more than a 
factor of 2 enhancement of ^'^S/^'^S in IRC+10216 while 
the other ratios are solar within error. An s-process en- 
hancement of ^'^Cl in IRC+10216 was found by Kahane 
et al. 1)2000(1 . These authors were able to match the obser- 
vation with the predictions of a low mass (1.5 Mq) AGB 
star model of solar metallicity. The enhancement predicted 
for the envelope after the 15th thermal pulse just cor- 
responded to the observed CI isotope ratio. Re-analysis, 
using the updated solar abundances, confirms the agree- 
ment, which persists for the last three thermal pulses with 
dredge-up (#15 to #17). Included in Fig. are the corre- 
sponding predictions for the sulfur isotopes (assuming, as 
in the case of CI, the updated solar starting composition). 
It is obvious in Fig.|31that qualitatively there is agreement 
between model and observations in that in IRC+10216 an 
enhancement of only '^^S/^^'S is seen. However, our ob- 
served enhancement of ~170% in this ratio is more than 
a factor of 4 higher than even the new prediction, which 
already is about twice as high (^40 % instead of ^20 %) 
than that in Kahane et al. (2000,) . 

It has been argued in the past that the central star of 
IRC+10216 might correspond to an AGB star of interme- 
diate mass, of around 5M0, instead of the low mass model 
discussed here. This possibility was at the upper limit of 
the range of uncertainty of the available distance estimates 
of 110 — 170 pc (Winters et al. IBMI Crosas & Menten 
[TM71 Le Bertrc 1997J Groenewegen et al. HMHl Weigelt 
et al. 1199 81. An intermediate mass was inferred by Guelin 
et al. |{T^95|) on the basis of older estimates of some iso- 
tope ratios (see the discussion in Kahane et al. 2000). We 
have computed a series of AGB models of different initial 
mass using updated network of neutron capture cross sec- 
tions and solar abundances. These confirm the conclusions 
already reached by Kahane et al. (2000) that an interme- 
diate mass of 5 Mq is excluded by the almost solar Mg iso- 
tope ratios observed. Indeed, in intermediate mass stars 
the ^^Ne(a,n)^^Mg and ^^Ne(Q;,7)^^Mg reactions would 
be more efficiently activated by the higher temperature 



in the convective He flashes, providing a clear excess of 
both neutron-rich Mg isotopes. Supporting evidence de- 
rives from the observed '^^Cl/'^^Cl ratio, which would be 
much higher in an intermediate mass star. The resulting 
34g^36g i-atio, in any case, does not depend much on the 
initial mass. 

While the enhancement of the observed ^^S/^^S ratio 
with respect to the 1.5 Mq models could, in principle, be 
due to some unrecognized problems with the AGB models 
or assumptions on nuclear parameters involved with it, we 
believe this to be highly unlikely and the main s-process 
in IRC+10216 to be not the main source of its (relative 
to solar) enhanced '^^S/'^'^S ratio. Alternative explanations 
are briefly discussed below. 

5.4.3. Chemical evolution 

One clue may be the close agreement between the value for 
IRC+10216 and the one reported for the ISM (also shown 
in Fig.EJ by Mauersberger et al. (1996), which at first 
glance is puzzling. A possible solution may be afforded by 
Galactic chemical evolution (GCE) and the different times 
when IRC+10216 and the Solar System formed, i.e. as- 
suming a starting composition different from solar. Based 
mostly on the ^^C/^^C ratio and limits to C/0, Kahane et 
al. { 2U^ derived an upper limit of 2 Mq for IRC+10216, 
which nevertheless still allows the object to be as young 
as ^1.5 Ga. Assuming a correspondingly evolved start- 
ing composition, observations might then be compatible 
(within errors) with model predictions. Indeed, first, there 
is no need for a large intrinsic contribution to ^^S from 
IRC+10216 (in agreement with the relatively small ef- 
fect predicted by the AGB star model; Fig.lSJ. And sec- 
ond, the large difference between the current ISM (and 
IRC+10216) and the Solar System, is compatible with 
massive star model predictions, where the '^^S Galactic 
contribution increases with metallicity (as already pointed 
out by Mauersberger et al.. 11996]) . This would require a 
rather different evolution of the CI and S isotopic ratios 
during the last few Ga, since the observed ^^Cl/^^'Cl ra- 
tio in IRC+10216 is matched by AGB model predictions 
using solar as the starting composition. Indeed, both CI 
isotopes are thought to be synthesized in massive stars in a 
primary way during explosive Ne and O burning (Weaver 
& Woosley 1995). Only 4% of the solar abundance of ^''Cl 
is contributed by the main s-component in low mass AGB 
stars, while '^^Cl gets actually depleted by neutron cap- 
tures, according to the GCE calculations by Travaglio et 
al. H2004|) . On the observational side, ^^Cl/^^Cl in Orion 
has been found to be enhanced relative to solar (Salez et 
al., I1996|l . but the errors are too large to reach a firm 
conclusion. 

A further possibility is that IRC+10216 simply started 
with a '^^S abundance that was unusually high for its time 
of formation and location, i.e. invoke inhomogeneity in the 
ISM. Of course, it could as well be the Solar System rather 
than IRC+10216 that started out with an unusual com- 
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position. Both suggestions are not purely ad hoc. Aheady 
Lugaro et al. (|1999J . in order to explain isotopic variations 
in the Si isotopes of single presolar silicon carbide grains 
from primitive meteorites, called for small local inhomo- 
geneities in the ISM at the time of parent star formation 
of the grains as an alternative to Galactic chemical evo- 
lution. And, as has been known for quite some time al- 
ready, oxygen isotopic ratios in the Solar System seem to 
be quite peculiar (e.g., Wilson & Rood, 119941 Prantzos 
et al.. ll99 ()ll. In addition, the interstellar Rb isotope ratio 
^^Rb/^'^Rb determined toward p Ophiuchi A (Federman 
et al. I2UI)4|I indicates a higher ratio of r- to s-process nu- 
clides in the Solar System as compared to the local ISM. 

Reifarth et al. H2(JU(J|I discussed various other possible 
nucleosynthesis scenarios besides the s-process. Accepting 
Galactic chemical evolution as the reason for the difference 
between the Solar System and the ISM ■^^S/^^S ratios, 
those that involve core collapse supernovae (SN) (other 
than the s-process) face the problem of how to account 
for the large increase during the last 4.6 Ga. This prob- 
ably includes the r-process. To some extent this also in- 
cludes the "neutron burst" , intermediate between s- and 
r-process and thought to occur during passage of a shock 
wave through the He shell of type II SN, which Meyer 
et al. H2000|l modelled to explain specific types of isotope 
abundance anomalies found in supernova grains preserved 
in meteorites. Note, however, that the efficiency of the pro- 
cess depends on previous s-process seeds and that its ef- 
fects are included in the SN II yields calculated by the au- 
thors cited above. Low-entropy scenarios associated with 
supernovae type la that can account for production of 
neutron-rich isotopes such as ^^Ca and ^"Ti (Meyer et al. 
119961 Kratz et al. 1200 l|l would be qualitatively compatible 
with the trend in GCE, but again existing information on 
production yields for ^^S is not encouraging (Reifarth et 
al.Hnnni Kratz et al. lMTTji . 

In summary, while it appears we have observationally 
identified a stellar source with enhanced '^^S and there are 
some reasonable ideas about how Galactic ^^S is produced, 
we are left with a number of open questions. Further stud- 
ies of other objects and/or regions in the ISM, especially 
when correlated with other isotopic systems such as the 
carbon, oxygen, and chlorine isotopes, may turn out to be 
helpful. The same holds true for continued study of the 
relevant stellar physics and nucleosynthesis processes. 

6. Conclusions 

The main conclusions of this paper are: 

1. Measuring the rare isotopic species of CS and SiS 
in IRC-hl0216, the abundance ratio of ^^S/^^?, in 
IRC-l-10216 is 107(±15), comparable to values in the 
Galactic interstellar medium (115), but smaller than 
the solar value (288). 

2. The increase of the '^^S abundance relative to '^^S only 
qualitatively follows the prediction of a low mass AGB 
star model. Quantitative agreement of the observed 



Table A.l. Line identifications 



Nr. 


Tl 


Freq." 


Im.-Freq" 


Identification 




mK 


MHz 


MHz 




1 


4 


94903 


98126 


HCCC"CN 36 — 35 


2 


3.5 


94913 


98155 


SiC4 32 - 31 (im.) 


3 


5.5 


95018 


98011 


I-C3H'' (im.) 


4 


6 


95035 


97095 


C^^S; I-C3H'' (im.) 


5 


21 


95049 


97980 


CS(2-1) (im.) 


6 


20 


95087 


97942 


SiC4 31-30 


7 


1.5 


U95115 


97914 




8 


1.5 


95126 


97903 


SiCCC^^C 32 - 31 


9 


20 


142297 


151263 


H2C4 16l,l6 — 152,15 


10 


23 


142321 


151235 


Al^'^Cl 10 - 9" 


11 


3.5 


142345 


151215 


HCCC"CN 54 - 53 


12 


5 


142402 


151158 


SiS J = 8 - 7,t; = 4"^? 


13 


29 


142410 


151150 


NaCN 927 - 826 


14 


5 


U142489 


151071 




15 


40 


142502 


151058 


CCS llii - 10^0 


16 


7 


142525 


151035 


C36s 2-2" 


17 


460 


142560 


151000 


^^SiS 8-7" 


18 


4 


U142594 


150966 




19 


12 


U142622 


150940 




20 


7 


142675 


150885 


SiCCC^C 48-47 


21 


9 


142695 


150865 


artefact 


22 


850 


142730 


150831 


C4H iV = 15 - 14= 


23 


9 


142755 


150805 


C4H iV = 15 - 14? 


24 


850 


142768 


150793 


C4H iV = 15 - 14^= 



a) The band from which the line probably originates is given in 
bold face; for information, the frequency of the corresponding 
image band is also given, b) Thaddeus et al. HM^ . c) 
Cernicharo et al. (2000) and references therein, d) The SiS 
« = 4J = 8- 7 line is at 142399.8 MHz. A shift of 1 - 2 MHz 
could be possible if the line is produced in the acceleration 
region (blue). 

^''S/^^S ratio with the stellar models can be reached 
if the age of IRC-f 10216 and Galactic Chemical evolu- 
tion are taken into account or if there are considerable 
inhomogeneities in the interstellar medium, and either 
IRC-l-10216 or the Sun started with an unusual ^^S 
abundance. 

3. From the observed line density toward IRC+10216 and 
toward Galactic star forming regions, confusion limits 
toward these sources are estimated. 

Appendix A: Line identifications 

In Table lXm we list the observed spectral line features 
seen in Fig.^ and the approximate values for TJ, the rest 
frequency and the rest frequency in the image band (as- 
suming an LSR velocity of — 26kms~^). The final column 
contains some identifications from the literature. In the 
256 MHz wide spectrum of the C^^S 2 — 1 line we have 
identified 8 spectral line features. Of those, three can be 
identified with transitions from the image sideband. In 
the 256 MHz around the frequency of the C^^S 3-2 
line we have detected 16 spectral line features, none of 
which can be assigned to a line from the image band. Six 
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Table B.l. The confusion limit reached by the IRAM 
30-m telescope toward some molecular line sources 



Source 
Av 


A 

mm 


weakest 
detect. Lines 
mK 


,a) 
"•int 

h 


Ref.'' 


IRC+10216 


3 


1.5 


28 


1) 


28kms"^ 


2 


3.5 


10 


1, 2) 




1.3 


^2 


80 


2) 


Orion-KL 


3 


20 


0.8 


3) 


5 kms~^ 


2 


50 


0.3 


4) 




1.3 


70 


0.4 


3) 


SgrB2(N) 
~ 17kms"^ 


3 
2 
1.3 


~ 20 
~ 30 
~ 40 


0.4 
0.4 
0.5 


3) 
3) 
3) 



a) integration time (on+off) with present receivers at the 
IRAM 30-m telescope under normal winter conditions (good 
summer conditions), one polarization to obtain an rms of 1/3 
Tmb with a velocity resolution of 1/5 the line width (one polar- 
ization receiver only); b)References: 1) this paper 2) Ziurys et 
al. (I2002t : 3) unpublished data; 4) Mauersberger et al. 119881 : 
c) weaker lines can be identihed if one makes use of the unique 
line shapes of the spectra of this source (Cernicharo et al. 

lines can be identified with known transitions from the 
Cernicharo et al. lj2()()0(l line survey, and the other iden- 
tifications were made with the catalog of the Instituto de 
la Estructura de la Materia (Madrid) containing ^ 1000 
species described by the authors. There is a clearly de- 
tected feature at the frequency of the C'^^S line; another 
feature at 142.755 GHz is tentatively assigned to emission 
from C4S. Eight further lines have not been identified. The 
emission at 142.695 GHz is an artefact due to some bad 
spectrometer channels. 

Appendix B: The confusion limit toward 

astronomical molecular line sources 

It is evident from Sect. 15. 1I that in our 2 and 3 mm spec- 
tra of IRC-l-10216 we are close to the confusion limit, 
where an increase of integration time does not yield much 
further information. It is interesting to investigate where 
line confusion begins to play a role for other favorites of 
molecule hunters such as the Orion Hot Core, SgrB2 or 
the starburst galaxy NGC253. We limit our discussion to 
the IRAM-30m telescope. If one wants to extrapolate our 
results to other telescopes with a higher or lower resolution 
one has to take into account the detailed source structure 
chemical and physical conditions within the regions ob- 
served (Comito & SchilkeEDUa. 

The definition of the useful observing time or rms to 
be reached is by no means unique and depends on whether 
one is interested in a mere detection of a line or whether 
one also wants to obtain some detailed information on 
the line shape. Here we try to be pragmatic: we have in- 
vestigated spectra of several molecular line sources made 
with the IRAM 30-m telescope at 3 mm, 2 mm and 1.3 mm 



wavelength made with such a long integration that ^ 50% 
of the spectral range observed is covered with features. We 
have determined the antenna temperature of the weakest 
unequivocal line features (knowing the typical line width 
and shape in the sources studied) . The results are given in 
Table l^TI We also give an indication of the necessary rms 
for a 3ct detection with a velocity resolution 1/5 the full 
line widths typical for these sources and the correspond- 
ing necessary integration time with the 30-m telescope and 
its present receivers (including all observing overheads, us- 
ing one polarization only) . In Orion, part of the confusion 
arises because of the large line widths observed in the out- 
fiow source. This can e.g. be prevented by observing at a 
carefully selected position offset from the outflow source 
(Combes et al. HM^ . 
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Fig. 2. C^^S and C^^S J = 2 - 1 and 3-2 spectra toward 
IRC+10216. The channel spacing corresponds to 3kms~^ 
for the 2 - 1 hnes of C^^S and C^^S, 2 km s'^ for the 3 - 2 
hne of C^^S and 4kms-i for the 3 - 2 transition of C^^S. 
Intensities are in mK (T^). Also shown are fits (convolved 
with the channel widths) of a shell type emission profile to 
the lines. The lineshapes of the C'^^S transitions were de- 
termined from the fits to the corresponding C'^^S lines. In 
the frame for C'^^S (2 — 1), the emission feature at the ex- 
treme left corresponds to the 2 — 1 line of the main isotope 
of CS observed in the image band. A fit to the C'^^S hne 
and to two I-C3II lines from the image band fixing the in- 
tensity ratios to the values obtained from a well calibrated 
spectrum (Mauersberger et al. I1996|l is overlayed on the 
spectrum. In the third frame the residuals obtained from 
subtracting the fit to the lines from the image sideband 
are shown. 
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Fig. 3. Si^'^S and Si^^S J = 5-4 and 6 - 5 spectra toward 
IRC-1-10216. The channel spacing corresponds to ~3.5 — 
4 kms^^. Intensities are in mK (T^). Also shown are fits 
(convolved with the channel widths) of a shell type emis- 
sion profile to the lines. In case of the Si'^^S J = 6 — 5 tran- 
sition, a simultaneous fit (assuming identical line width) 
with a blended unidentified line is shown. 
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Fig. 4. Isotopic ratios of '^^S/'^*S obtained from four 
transitions. The straight line indicates the weighted mean 
value; a dashed line indicates the error of the mean value. 
Also shown is the solar system ratio (3.47 10""^). 
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Fig. 5. Sulfur isotopic ratios in IRC+10216, normalized to 
■^"'S and solar ratios. Data from Kahanc ct al. ( 2000J and 



this work. Also shown are the S isotope ratios predicted 
for the envelope of a 1.5 M© solar-metallicity, TP-AGB 
star at the time the CI isotope ratio is matched (see text; 
with updated solar abundances), and the ratios measured 
for the ISM by Mauersberger et al. 



